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ABSTRACT : The characterization of the electromagnetic interaction for 
a spacecraft in the lunar environment, and identification of viable charging 
mitigation strategies, is a critical lunar mission design task, as spacecraft 
charging has important implications both for science applications and for 
astronaut safety. To that end, we have performed surface charging 
calculations of a candidate lunar spacecraft for lunar orbiting and lunar 
landing missions. We construct a model of the spacecraft with candidate 
materials having appropriate electrical properties using Object Toolkit and 
perform the spacecraft charging analysis using Nascap-2k, the 
NASA/AFRL sponsored spacecraft charging analysis tool. We use 
nominal and atypical lunar environments appropriate for lunar orbiting and 
lunar landing missions to establish current collection of lunar ions and 
electrons. In addition, we include a geostationary orbit case to 
demonstrate a bounding example of extreme (negative) charging of a lunar 
spacecraft in the geostationary orbit environment. Results from the 
charging analysis demonstrate that minimal differential potentials (and 
resulting threat of electrostatic discharge) occur when the spacecraft is 
constructed entirely of conducting materials, as expected. We compare 
charging results to data taken during previous lunar orbiting or lunar flyby 
spacecraft missions. 


1- INTRODUCTION 

Development of space systems for reliabile operation in lunar environments will necessarily need to 
consider spacecraft charging. It has been known since the first exploration of the Moon by robotic 
and manned spacecraft that the lunar dayside charges to a few tens of volts where the photoelectron 
current is the dominant charging process [1,2] and negative nightside potentials develop due to the 
preferential collection of hot electrons that penetrate the plasma wake [3]. The early observations 
and theoretical studies are consistent with the more recent results from the Lunar Prospector 
spacecraft which demonstrate positive lunar surface potentials in the daylit hemisphere and negative 
surface potentials reported on the order of a few hundred volts in the lunar wake [4, 5, 6, 7] or even 
extreme values exceeding a few kilovolts during periods when solar disturbances provide an 
additional source of hot electrons [8], Finally, we note that the analyses presented in the Lunar 


Prospector papers confine their results to the potential of the lunar surface and no attempt is made to 
evaluate the charging properties of the spacecraft itself in orbit about the Moon. 

Spacecraft charging in lunar environments is not fundamentally different than conditions 
encountered in previous missions within the Earth’s magnetosphere and in the near Earth 
interplanetary space. The renewed interest in lunar exploration and a return of humans to the Moon 
however places a particular emphasis on designing safe and reliable systems for operations under a 
variety of charging conditioas. Many of the spacecraft charging tool sets currently in use for 
evaluating potential distributions and electric fields on the surface and in the space surrounding a 
spacecraft focus on low Earth orbit (including auroral charging), geostationary orbit, and 
interplanetary space where the majority of space vehicles are located. It is important to understand 
if these tool sets are adequate to describe the low density, high temperature environments of the 
lunar wake as well as the variety of conditions experienced by a spacecraft in lunar orbit while the 
Moon transits the Earth’s magnetotail. 

This paper presents Nascap (NASA and Air Force Charging Analyzer Program) -2k [9] surface 
charging analyses for a candidate lunar orbiting spacecraft to determine if the current version of the 
Nascap-2k 3-D charging analysis model is adequate for use in designing unmanned lunar orbiting 
spacecraft and the transportation systems required for human missions to the Moon. The current 
default Nascap-2k options for computing spacecraft potentials include interplanetary, geostationary, 
and low Earth orbit environments. Since there is no lunar specific option, we use both the 
interplanetary and geostationary orbit options here for computing spacecraft surface potentials and 
electric fields for a range of environments from solar wind, magnetosheath, and into the Earth’s 
magnetotail. Discussion of the results include an assesment of which of the existing options in 
Nascap-2k provide the most applicable results for lunar applications and what modifications would 
be desirable to support the analyses required to design space systems for upcoming NASA missiosn 
to the Moon. 


2 - MODEL 

A sample model for the Lunar Prospector was written in Object Tool Kit (OTK), the object editing 
module of Nascap-2k. Figure lb shows the completed OTK Lunar Prospector model based on an 
artist rendition of the vehicle given in Figure la. Care was taken to represent as acuarately as 
possible the actual Lunar Prospector spacecraft geometry, however, there is no expectation given 
that all materials and dimensioas are correct. The body of the Nascap- 2k spacecraft (cylinder) is 
covered in solar cells. The top and bottom of the body is covered in graphite. There are three 
booms extending radially outward from the body. For Nascap-2k purposes, they have been 
modeled as graphite. At the end of each boom are the science modules. Referring to Figure 1 b, the 
alpha particle spectrometer is represented as the square module in the lower right hand comer of the 
image. It is modeled as graphite with aluminum square patched on all sides except the side 
connecting to the boom. The cylinders extending outward from the alpha particle spectrometer are 
the neutron spectrometers: one is shown as tin and the other as cadmium. The gamma ray 
spectrometer is extending back on the right side of the image. It is modeled as graphite. The 
electron reflectometer is attached to the boom extending to the back and left in Figure lb. The 
magnetometer is then attached to it via another boom. Both the electron reflectometer and the 
magnetometer are modeled as graphite. For simplicity, tin and cadmium were created with the 
same default Nascap-2k properties as aluminum. All remaining materials also use the Nascap-2k 
default properties. 



a. b. 

Figure 1 - (a) Artist’s rendition of the Lunar Prospector spacecraft showing spacecraft bus covered with 
solar arrays and booms for instrumentation (image courtesy http://nssdc.gsfc.nasa.gov/planetarv/lunarprosp.hunl) . 
(b) Nascap-2k Object Tool Kit model of the Lunar Prospector spacecraft with surface materials adopted for 
charging analysis indicated by the color scheme including graphite (grey), solar cells (dark blue), cadmium 

(purple), tin (yellow) and aluminum (light blue). 


3 - ENVIRONMENT 

The Moon spends approximately 25% of the time inside the Earth’s magnetosheath and 
magnetotail, with the rest of the time spent inside the solar wind. The Moon has no appreciable 
atmosphere and no global magnetic field, so the charged particles comprising the solar wind, 
magnetosheath, and magnetosphere interact directly with the lunar surface. Plasma flows produce a 
plasma wake on the downstream side of the Moon where density depletions within the wake of 
some two to three orders of magnitude are observed [7], The first observations of the lunar wake 
were made by Explorer 35 and the Apollo sub-satellites [10]. There have been a number of more 
recent observations with spacecraft such as WIND [11] during lunar fly-by maneuvers and Lunar 
Prospector [7] in low lunar orbit. 


In order to estimate charging conditions within the wake, plasma environment parameters relevant 
for the wake conditions are required for input to the Nascap-2k charging model. We use an 
analytical model of the wake plasma density and temperature environments for specifying the 
charging environments in the lunar wake. The wake model is based on functions describing plasma 
density and temperature variations as a function of depth into the lunar wake and distance from the 
Moon derived by Halekas et al. [7] from Lunar Prospector electron temperature and density 
observations and model ion environments derived by Samir et al. [ 12]. Based on the ambient solar 
wind or magnetotail plasma environment, the model calculates the disturbed density and 
temperature in the wake region for electrons, protons, and alpha particles. 


Free-field (upstream) input parameters for the wake model were taken from Paterson and Frank [13] 
and Feldman et al. [14] for the magnetotail boundary layers, Feldman et al. [14] for the 
magetosheath. Craven [15] and Feldman et al. [14] for the plasma sheet, and from Newbury et al. 
[16] and Feldman et al. [14] for the high and low speed solar wind environments. The wake model 
then provides the plasma conditions at 100 km altitude in the deep lunar wake. Environments given 
in Table 1 represent the model results for the deep lunar wake which are used as the Nascap-2k 
input parameters for both interplanetary and geosynchronous (GEO) runs. Slight modifications in 


the lunar environment model output may have been made to some, but not all parameters to 
correspond to the inputs needed for Nascap-2k. 

Table 1. Environments 



Boundary Layer 

Magnetosheath 

Plasma Sheet 

Solar Wind High Speed 

Solar Wind Low Speed 


Inter- 


Inter- 


Inter- 


Inter- 


Inter- 



planetary 

Geo 

planetary 

Geo 

planetary 

Geo 

planetary 

Geo 

planetary 

Geo 

P(m 3 ) 

7.326x10 4 


6.64x10 5 


2.428x1 0 4 


266.3 


6374 


T.(eV) 

156.7 

156.7 

290.9 

290.9 

1316 

1316 

78.28 

78.28 

47.78 

47.78 

Vi (m/s) 

5.4x10* 


4x1 0 a 


2x10° 


7.02x10" 


3.27x10° 


& (eV) 
P« (m^ 1 ) 

15.22 

7.326x10* 

835.3 

6.64x1 0* 

208.8 

2.428x10* 

2573 

266.3 

558.2 

6374 

P^rn 3 ) 


6.684x10* 


6.06x10" 


2.215x10* 


50 


50 

T,(eV) 


111 6 


920 


995.7 


10 


10 

U (A/m 2 ) 

2458x10* 

2.458x10* 

3.036x10' 

3.036X10* 

2 361x10* 

2.361x10* 

6.316x10 n 

6.316x10 n 

1.181x10* 

1 181x10* 

MA/m 2 ) 

6.338x10 10 

4416x10 10 

4.255x10* 

1.15x10* 

7.779x10 10 

4.372x10* 

2.995x10 11 

9.891x10 14 

3.339x10 10 

9 891x10 14 

Charging time 

300 sec 

300 sec 

1000 sec 

300 sec 

1000 sec 

300 sec 

1000 sec 

300 sec 

1000 sec 

300 sec 


4 - RESULTS 

The following subsections describe the results of the Nascap-2k interplanetary and geosynchronous 
runs for boundary layer, magnetosphere, plasma sheet, high speed solar wind, and low speed solar 
wind environments. The problem type is a surface charging analysis using analytic currents. All 
runs were for at least 300 seconds of charging time, with no grid, no magnetic field or sun, an initial 
electron particle species, and an initial applied potential of -5 volts. 


4.1 - BOUNDARY LAYER 


Figure 2 shows the potential results of the interplanetary and geosynchronous runs for the boundary 
layer calculations. Potential results were exactly the same in both runs : 1 .60 to 2.90 volts, with the 
solar cells at 2.90 volts and the graphite areas of the spacecraft at 1 .6 volts. For the potential results 
the interplanetary and geosynchronous runs went to equilibrium immediately and stayed smooth. 
However, the charging current and electric field results had numerical noise for the geosynchronous 
run. Interestingly, when plotting the electric field for the interplanetary run, results stabilized at 
zero V/m. However, for the geosynchronous run, the electric field was between -350 and 400 V/m. 



a. 


b. 


Figure 2 - Nascap-2k boundary layer environment potential results using the (a) interplanetary and (b) 

geosynchronous potential computation options. 



4.2 - MAGNETOSHEATH 


Figure 3 shows the potential results of the interplanetary and geosynchronous runs for the 
magnetosheath. Potential results were slightly different for each run : 1.60 to 4.80 volts for the 
interplanetary run and 1.0 to 3.4 volts for the GEO run. The solar cells showed more variation in 
end results in the interplanetary run than the GEO run. Again, graphite charged lower than the solar 
cells for both cases. For the potential results the interplanetary and geosynchronous runs went to 
equilibrium immediately and stayed smooth. However, the charging current and electric field 
results for both had small numerical noise. When plotting the electric field for the interplanetary 
run, results stabilized at zero V/m. However, for the geosynchronous run, the electric field was 
between -150 and 400 V/m. 



a. b. 

Figure 3 - Nascap-2k magnetosheath environment potential results using the (a) interplanetary and (b) 

geosynchronous computation options. 


4.3 - PLASMASHEET 

Figure 4 shows the potential results of the interplanetary and geosynchronous runs for the 
plasmasheet calculations. Potential results were -2.5 to 3.5 V for the interplanetary run and -1.5 to 
3.0 V for the GEO run. Again, the solar cells for the interplanetary run showed more variation in 
charging. Graphite again was more negative. For the potential results the interplanetary run went 
to equilibrium immediately and stayed smooth, while the GEO run had a sizable amount of 
numberical noise around equilibrium. There was a small amount of noise for the charging current 
for both types of runs. The electric field for the GEO run was very noisy around the minimum and 
maximum values. When plotting the electric field for the interplanetary run, results stabilized at 
zero V/m. However, for the geosynchronous run, the electric field was between -800 and 800 V/m. 




Figure 4 - Nascap-2k plasmasheet environment potential results using the (a) interplanetary and (b) 

geosynchronous computation option. 


4.4 - SOLAR WIND - HIGH SPEED 

Figure 5 shows the potential results of the interplanetary and geosynchronous runs for the high 
speed solar wind calculations. Potential results for the interplanetary and geosynchronous runs 
were almost identical at 0.78 to 0.81 V for interplanetary and 0.79 to 0.80 V for geo. Both went to 
equilibrium almost immediately and were smooth. The solar cells charged more positive than the 
graphite for both cases. The charging current and electric field were the same for both with the 
electric field between -200 and 250 V/m. 



a. b. 

Figure 5 - Nascap-2k high speed solar wind environment potential results using the (a) interplanetary and (b) 

geosynchronous computation option. 


4.5 - SOLAR WIND - LOW SPEED 

The low speed solar wind case had the largest negative charging of any of the cases ran. Figure 6 
shows results from the interplanetary and GEO runs. However, there was basically no differential 
charging. The interplanetary run produced a potential of -235.6 V, while the GEO run was -154 to 
-153.9 V. The potential, charging current, and electric field profiles for both interplanetary and 






GEO runs were smooth and went to equilibrium within the final charging time. The electric field 
for the interplanetary case was -8xl0 4 to 6xl0 4 V/m. The eletric field (E-field) for the GEO case 
was -5xl0 4 to 4xl0 4 V/m. 



Figure 6 - Nascap-2k low speed solar wind environment potential results using the (a) interplanetary and (b) 

geosynchronous computation option. 


Nascap-2k results for all the environments for both interplanetary and geosynchronous runs are 
summarized in Table 2. 

Table 2. Charging Results 



Boundary Layer 

Magnetosheath 

Plasma Sheet 

Solar Wind High Speed 

Solar Wind Low Speed 


Inter- 

planetary 

Geo 

Inter- 

planetary 

Geo 

Inter- 

planetary 

Geo 

Inter- 

planetary 

Geo 

Inter- 

planetary 

Geo 

0>(V) 

1 60 to 2.90 

1.6 to 2.9 

1.6 to 4 80 

1 to 3.4 

-2.5 to 3 5 

-1.5 to 3 

0 78 to 0.81 

0.79 to 0 8 

-235.6 

-154 to - 
153.9 

equil 

immed 

noise -2V 

immed 

small noise 

immed 

very noisey 

immed 

immed 

to equil 

to equil 

E (V/m) 

-0 

-350 to 400 

-0 

-150 to 400 

-0 

-800 to 800 

-200 to 250 

-200 to 250 

-8x104 to 
6x104 

-5x104 to 
4x104 

equil 

Charging 

-1.2x10-6 to 

noise 

-1.2x10-6 to 

-12x10-6 to 

small noise 
-1.0x10-6 to 

-12x10-6 to 

min and 
max very 
noisey 

-12x10-6 to 

immed 
-1.1x10-6 to 

immed 
-1.1x10-6 to 

to equil 
-1.1x10-6 

to equil 
-11x10-6 to 

Current (A/m*) 

1.2x10-6 

1.2x10-6 

1.2x10*6 

12x10-6 

12x10-6 

1 2x10-6 

1.1x10-6 

1.1x10-6 

to 1.1x10-6| 

| 1 1x10-6 

equil 

Charging time 

300 sec 

noise -0 
300 sec 

small noise 
-O 

1000 sec 

tiny noise -0 
300 sec 

1000 sec 

tiny noise 
300 sec 

mmed 
1000 sec 

immed 
300 sec 

to equil 
1000 sec 

to equil 
300 sec 


5 - CONCLUSION 

Differences in potentials which result from using the interplanetary and geosynchronous orbit 
computational options are likely due to how Nascap-2k treats the plasma current for interplanetary 
problem types as opposed to geosynchronous problem types [Mandell et al., ] since the 
computational techniques are similar for the both, due to both having a long Debye length and 
tenuous plasma, except for the plasma current treatment. 





While the rule of thumb for charging in darkness is that a vehicle will reach an equilibrium potential 
of a few kT, the results obtained for the lunar wake suggest the majority of the conditions result in 
positive potentials. It is possible to reconcile the positive (or small negative) spacecraft potentials 
obtained in the Nascap-2k charging results with the 0.2 to 1.4 kV negative potentials reported by 
Halekas et al. [6, 7, 8] for the lunar surface by considering that the default secondary electron yields 
incorporated in Nascap-2k for the solar cell material approach six for kilovolt electrons. The 
moderate electron temperatures used in the analyses reported here yield electron spectra dominated 
by electron energies of a few keV or less since Nascap-2k utilizes Maxwellian distribution functions 
for computing electron flux in both the geostationary and interplanetary analyses models. As such, 
the strong secondary yields of the solar cells covering a large surface area of the candidate 
spacecraft dominate the charging process and positive potentials result even in the lunar wake 
environments. Indeed, it is well known that temperature threshold effects exist for the onset of 
negative vehicle charging due to energy at which the secondary electron yield curve decreases 
below unity [19, 20]. Had more severe environments with larger electron temperatures been 
included in the case studies, there would have been more negative potentials resulting from the 
Nascap-2k charging results. 


It is likely the environments measured by the Lunar Prospector spacecraft were better represented 
by Kappa distribution functions [7] since the electron environments in the deep lunar wake are 
characterized by energetic electrons preferrentially filling the plasma void ahead of the cold 
electrons and ions. It is recommended to include options for future versions of the Nascap-2k 
software to model current collection processes with Kappa functions to evaluate this effect, a 
feature not available in the current version of the code. 


Programs such as NASA’s new Constellation lunar program would benefit from modification of the 
current version of the Nascap-2k surface charging model to include a dedicated lunar charging 
module which includes non-thermal electron and ion distributions, plasma flow velocities 
independent of the solar illumination direction, and other effects unique to the lunar environment. 
However, for spacecraft orbiting the moon, the interplanetary environment charging run does an 
adequate job in calculating the potentials. The GEO problem type seems to do a better job 
calculating the electric field however. But the GEO problem type has more numerical noise that 
may or may not be an issue. 
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Lunar Wake Model 
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Charging Environments for Nascap-2k Analysis 
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Magnetosheath Environment Results 
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Interptanetary module results Geosynchronous module results 



Interplanetary module results Geosynchronous module results 




Interplanetary module results Geosynchronous module results 


Interplanetary Results 


Interplanetary Results 


Interplanetary Results 

Spacecraft Solar Arrays*.: 4 236 V 


Spacecraft Solar Arrays*,: 2.109 V 


Spacecraft Solar Arrays,.: -235.6 V 

Booms: 1.762V 


Booms: -2478 V 


Booms: -235.6 V 

Alpha par tele detector 1.762 V 


Alpha per tele detector -2478 V 


Alpha par tele detector -236.6 V 

Neutron Spectrometer i 762 V 


Neutron Spectrometer: -2478 V 


Neutron Spectrometer: -235.6 V 

Gamma Ray Spectrometer 1.762 V 


Gamma Ftey Spectrometer -2478 V 


Gamma Ray Spectrometer 236 6 V 

Electron Reftectometer: 1 762 V 


Election Reftectometer -2478 V 


Electron Reftectometer -235.6 V 

Magnetometer: 1.762 V 


Magnetometer: 2478 V 


Magnetometer. -236.6 V 

Sunlght..: 1 802 V 


Suntght,. 2405 V 


Suntght..: -236.6 V 

Eclpse... 3.110 


Eclpse,.: 0.2163 V 


Eckpse,.: -235 6 V 




Summary of Surface Potentials 











